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ABSTRACT 

We present radial velocity (RV) measurements of our sample of 40 M dwarfs from our planet search programme with VUT+UVES 
begun in 2000. Although with our RV precision down to 2-2.5 m/s and timebase line of up to 7 years, we are capable of finding planets 
of a few Earth masses in the close-in habitable zones of M dwarfs, there is no detection of a planetary companion. To demonstrate this 
we present mass detection limits allowing us to exclude Jupiter-mass planets up to 1 AU for most of our sample stars. We identified 
6 M dwarfs that host a brown dwarf or low-mass stellar companion. With the exception of these, all other sample stars show low RV 
variability with an rms < 20 m/s. Some high proper motion stars exhibit a linear RV trend consistent with their secular acceleration. 
Furthermore, we examine our data sets for a possible correlation between RVs and stellar activity as seen in variations of the Ha line 
strength. For Barnard's star we found a significant anticorrelation, but most of the sample stars do not show such a correlation. 
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1. Introduction 

High-precision differential radial velocity (RV) measurements of 
stellar reflex motions induced by an orbiting companion have so 
far been the most successful method to discover extrasolar plan- 
ets and to characterise their orbital properties. Originally, RV 
planet search programmes have largely concentrated on main 
sequence stars of spectral types late-F through K. That only a 
comparatively small number of M dwarfs were included comes 
from their faintness, which requires large telescopes to perform 
high-precision RV measurements of a few m/s. For an under- 
standing of the formation and abundance of extrasolar planets it 
is important to determine the presence and orbital characteristics 
of planets around stars of as many different types as possible, and 
especially around this most abundant type of star 

Even if M dwarfs are faint and require large telescopes, how- 
ever, they have two advantageous characteristics when search- 
ing for terrestrial exoplanets in the habitable zones (HZ) with 
radial velocity (RV) methods. Compared to solar-like stars, (i) 
they have a lower mass (M = 0.1 - O.6M0), and (ii) the habitable 
zone is close-in aroun d this cooler and less luminous type of star 
(L = 0.0008 -O.O6L0,lGuinaniLEngle"2OO9). By habitable zone 
we understand the reg ion that allow s liquid water on the planet 
surface as described in iKasting et al.l(ll993l) . 



Table 1. M dwarfs with known planets and their masses discov- 
ered with the RV method. 



Star 


b 


c d 


e 


References 


GJ 876 


2.53Mj„p 


0.79Mj„p 7.53Me 




[1][2][3][4] 


GJ581 


15.6Me 


5.1Me 8.2Ma, 


1.9Me 


[5] [6] [7] 


GJ436 


21Me 






[8] 


GJ317 


0.71Mj„p 






[9] 


GJ 674 


ll.lMa, 






[10] 


GJ 849 


0.82Mj,p 






[11] 


GJ 176 


8.4Me 






[12][13] 


GJ 832 


0.64Mj„p 






[14] 



Note: The masses are minimum masses with the exception for GJ 436 
an d GJ 876. [1 ] Delfosse et al. 1998, [2,3] Marcv et al. 1998, 200ll 
[4] iRivera et al. 2005, [5] Bonfils et al. 2005, [6] Udrv et al. 2007, [7] 
Mavor et al. 2009, [8] Butler et al. 2004, [9] Johnson et al. 2007, [10] 
Bonfils et al. 20 07,, [1 1] Sutler et aL ,2q06i. [ 12] iForveille et al.„2009i. 
[13] Butler et all2009l [14] Bailev et alj2009l . 



For M dwarfs, the HZ is typically 0.03 - 0.4 AU. The RV 
amplitude induced due to a planet by the Doppler effect is 



K 



m sin I 



1 - yl{W+m)a 



= 28.4m/s.^(^^r' 



Ml 



Jup 



M a I 



* Based on observations collected at the European Southern 
Observatory, Paranal Chile, ESO programmes 65.L-0428, 66.C-0446, 
267.C-5700, 68.C-0415, 69.C-0722, 70.C-0044, 71.C-0498, 072.C- 
0495, 173.C-0606, 078.C-0829. 

** Radial velocity data are available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/| 



(1) 

It increases for closer distances a (shorter periods) and lower 
stellar mass M. Thus the RV amplitude induced by a planet in 
the HZ of an M dwarf is higher than that of a solar-like star 
M dwarfs are also ideal targets for astrometric follow-up due to 
their lower mass, as well as for transit observations. In combina- 
tion with the RV method, astrometry allows the resolution of the 
sin /-ambiguity and true masses to be obtained. 

At present there are only a few M dwarfs known to have plan- 
ets, as summarized in Table [T] The M4V star GJ 876 has two 
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Jovian planets orbiting in a 2:1 resonance dMarcv etal.ll200U 
see also lBenedict et al]|2002l for an astrometric determination of 
the mass of the ou ter planet), and a th ird planet in this system 
has been found by iRivera et al] (l2005h . The M3V star GJ 581 
is another multiple system with three Neptune-type planets. The 
latest d etection of a fourth planet in this system bv lMavor et al.l 
(12009^ is a new highlight in planet search with RVs. With only 
1.9Me, GJ 581 e is the lowest-mass planet found so far with 
the RV method. The same work also resulted in a revised pe- 
riod for the outer planet GJ 581 d, placing it inside the habitable 
zone. The planet arou nd the M2.5 V star GJ 436 first discovered 
with the RV method (Butler et al. 2004) also turned out to be a 
transiting one tGillon et al. 2007). So far it is the only known 
transit ing planet around an M dwarf. Most recently. iBailev et aU 
(|2009|) has discovered the first long-period planet around an 
M dwarf (GJ 832). A planet with a mass of 24.5M^ (P = 
10.24 d) around GJ 176 announced by Endl et al. (2008) is re- 
jected bv lForveille et alj (12009.) and Butler et al. (2009.) . Instead, 
both groups find evidence of another plane t with a shorter period 
(8.78 d) a nd a minimum mas s of 8.4Me dForveille et alj|2009l) 
and 12M„ lButler et al.l ( l2009h . 

Several radial velocity (RV) surveys of M dwarfs have re- 
sulted in few or no detections, indicati ng a lower f requency of 
planets compared to solar-like stars; e.g.. lEndl et al.l (2006) mon- 
itored 90 M dwarfs (including the first data for 21 stars from 
our sample) wit hout a planet detection. The sample studied by 
ICumming et al.l (^008) surveying 1 10 M dwarfs contained only 
two planet hosting M dwarfs (GJ 876 and GJ 436). HARPS 
guaranteed time project #3 had 50 nights on 120 M dwarfs 
within 1 1 pc, and so far it has revealed only three planet-hosting 
M dwarfs (GJ 876, GJ 581, and GJ 674). 

M dwarfs are an ideal testing ground for competing models 
of the formation of gas giants. While the classical core-accretion 
model has severe problems with forming Jupiter-mass planets 
in th e less massive p rotoplanetary disks even around M dwarfs 
(e.g. iLaughlin et al.l 12004) . the competing gravitational insta- 
bility model can also efficieri tl y form Jovian-t ype companions 
around M dwarfs (lBossll2006h . llda & Liiil (^2005) even predict a 
higher frequency of icy giant planets with masses comparable to 
Neptune in short periodic orbits for M dwarfs than for G type 
stars. 

Recent results of microlensing s urveys (e. g. OGLE-2005- 
BLG-390b, 5.5 M^. iBeaulieu et alJ 12006: OGLE-2005-BLG- 
169b. lGouldetani20()6r presumed M dwarfs) may also indicate 
that low-mass planetary companions might be abundant around 
M dwarfs. 



2. Targets and Observations 

Our sample consists of 40 M dwarfs and one M gianQ, which are 
listed in Table |2l All M dwarfs are brighter than V < 12.2 mag 
and nearby within a distance of 37 pc (33 M dwarfs even within 
20 pc). Their spectral types range from MO to M5. The stellar 
masses M were deriv ed from the mass-luminosity relation by 
iDelfosse et al.l (12000') using the absolute brightness in K band. 
As an indicator of activity we selected the X-r ay luminosity 
as d etected in the ROSAT all-sky survey (from iHiinsch et alj 
Il999>) . Detected as X-ray sources are Barnard's star (L^ - 0.1), 
GJ 1 (0.6), GJ 190 (0.9), GJ 229 (1.3), and the most active 



' GJ 4106 is listed in the Catalogue of Nearby Stars with a parallax 
of84±17mas(11.9pc). However, with the Hipparcos parallax G J 4 1 06 
should be an M giant. 



Table 2. Targets with their sp ectral type, visual magnitude V, 
distance d dvan Leeuwenl2007l) . and stellar mass M derived from 
the K-band mass-luminosity relation by Delfosse et al. (2000). 



Star 


Spec Type 


V [mag] 


d[pc] 


M[Mo] 


Barnard (GJ 699) 


M4Ve 


9.54 


1.82 


0.16 


GJ 1 


Ml. 5 


8.57 


4.34 


0.45 


GJ27.1 


MO. 5 


11.42 


23.99 


0.53 


GJ 118 


M2.5 


10.70 


11.65 


0.36 


GJ 160.2 


MOV 


9.69 


23.12 


0.69 


GJ 173 


Ml. 5 


10.35 


11.10 


0.48 


GJ 180 


M2 


12.50 


12.12 


0.43 


GJ 190 


M3.5 


10.31 


9.27 


0.44 


GJ218 


M1.5 


10.72 


15.03 


0.50 


GJ229 


M1/M2V 


8.14 


5.75 


0.58 


GJ263 


M3.5 


11.29 


16.02 


0.55 


GJ357 


M2.5V 


10.85 


9.02 


0.37 


GJ377 


M3 


11.44 


16.29 


0.52 


GJ422 


M3.5 


11.66 


12.67 


0.35 


GJ433 


Ml. 5 


9.79 


8.88 


0.48 


GJ477 


Ml 


11.08 


18.99 


0.54 


GJ510 


Ml 


11.05 


16.74 


0.49 


GJ620 


MO 


10.25 


16.44 


0.61 


GJ637 


M0.5 


11.36 


15.88 


0.41 


GJ682 


M3.5 


10.96 


5.08 


0.27 


GJ739 


M2 


11.14 


14.09 


0.45 


GJ817 


Ml 


11.48 


19.17 


0.43 


GJ821 


Ml 


10.87 


12.17 


0.44 


GJ842 


MO. 5 


9.74 


11.99 


0.58 


GJ855 


M0.5 


10.74 


19.15 


0.60 


GJ891 


M2V 


12.20 


16.08 


0.35 


GJ911 


MOV 


10.88 


24.26 


0.63 


GJ 1009 


Ml. 5 


11.16 


17.98 


0.56 


GJ 1046 


M2.5+V 


11.62 


14.07 


0.40 


GJ 1100 


MO 


11.48 


28.93 


0.57 


GJ 3020 


M2.5 


11.54 


22.78 


0.62 


GJ 3082 


MO 


11.10 


16.56 


0.47 


GJ 3098 


M1.5Vk: 


11.21 


17.86 


0.50 


GJ 3671 


MO 


11.20 


17.74 


0.50 


GJ 3759 


MIV 


10.95 


16.97 


0.49 


GJ 3916 


M2.5V 


11.25 


15.10 


0.49 


GJ 3973 


M1.5Vk: 


10.94 


18.23 


0.54 


GJ4106 


M2 


10.82 


110.50 


0.55 


GJ 4293 


M0.5 


10.90 


25.06 


0.57 


HG 7-15 


MIV 


10.85 


37.31 


0.78 


ProxCen (GJ551) 


M5.5Ve 


11.05 


1.30 


0.12 



Proxima Cen (1.7). The rest of the sample was not detected by 
ROSAT, implying that these stars are inactive. 

The observations started in 2000, initially with 20 targets. 
Typically, three consecutive spectra per night were taken for 
each star with exposure times of 90 - 900 s depending on the 
object brightness. In April 2004 (JD=2453100) the sample was 
complemented by 21 additional stars, while the monitoring of 
HG 7-15 was ended. Since then the number of spectra per 
night was reduced to one, with the exception of Barnard's star, 
GJ 160.2, GJ 821, and Proxima Cen. 

Our data of the first 20 stars (-1- GJ 510) taken before mid- 
2005 were already included in the study by Endlet al. (200^. 
Here we present the full data set for all 41 stars as observed 
until March 2007. A d etailed study of the fu ll Proxima Cen 
data set can be found in lEndl & Kiirsteil (l2008l) . A study of the 
pre-2002.75 data set of Barnard's star is given by iKiirster et al.l 
(2003). 

The observations were carried out with the UVES spectro- 
graph at the VLT-UT2, directly fed via image slicer #3 that re- 
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Fig. 1. Radial velocities for 35 M dwarfs. Dashed lines represent no fit; they show the predicted secular acceleration effect caused 
by the proper motion (only for Proxima Cen and stars with > 1 m/s/yr ). 
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Fig. 2. Radial velocities for 6 more M dwarfs. The orbital solution for GJ 1046 is taken from lKiirster et alJ (l2008l) . 



distributes the light from a l"xl" aperture along a 0.3" wide 
slit. This resulted in a resolving power of R - 100 000 - 120 000. 
The red arm of UVES was employed and a wavelength coverage 
495 - 704 nm over 37 orders was selected. An iodine cell was 
used for precise wavelength calibration and modelling of the in- 
strumental profile. Only the range of 500 - 600 nm was used to 
determine the RVs. This range is rich in iodine lines. 

The data reduction included bias subtraction, flat-fielding, 
Echelle straylight subtraction , wavelength calibra tion, and 
barycentric correction (see also lEndl & Kurstej|2008l) . The data 
modelling with the "A USTRAL" code to obtain the RV is de- 
scribed in lEndl et alj (12000). We achieve an RV precision of 
2m/s for bright stars. In practice, photon noise limits the RV 
precision for faint M dwarfs and correspondingly the errors are 
larger for those stars that were obse rved with lower S/N . Radial 
velocities and also their errors (see iKiirster et al.ll2003l for dis- 
cussion) were combined into nightly averages. They are shown 
in Figs.[T]and|2](the RV data are available as online material). 



3. Data analysis 

3.1. Secular acceleration 

Even if a star moves undisturbed with a constant space velocity 
V, it can show a change in its radial velocity (RV) . This secular 
or per spective acceleration was first measured by IKiirster et al.l 
( l200 3l) for Barnard's star In polar coordinates the radial and tan- 
gential velocity component, respectively, are (Fig.[3]l 

Vr — —v cos tp and Vt = vsini;o. 

The diff'erentiation of Vi with respect to time f yields the secular 
acceleration (v = const.) 

. dv,(f) 

Vr = - V Sm (f ■ (fi - Vt(p. 




observer 



Fig. 3. Change in the radial velocity Vr due to constant motion v 
(secular acceleration v^). 

The tangential velocity Vt - t/ // depends on the distance d of the 
star and its proper motion which can also be identified with 
the time derivative of the angular position = //. Therefore the 
instantaneous secular acceleration is given by 



dvAt) 
dt 



7=^ 



-d = 22.98 



m/s (ju^ + /"a) ■ yrVai-csec^ 



yr 



;r/mas 



(2) 



where //„ and /us are the proper motion in right ascension and 
declination, respectively. It only depends on the proper motion fi 
and parallax n, which are easily accessible from the Hipparcos 
catalogue. Note that the knowledge of Vi and v is not explicitly 
required for the prediction of the instantaneous secular accelera- 
tion. 

One has to take this effect into account, especially for high 
proper motion stars, to avoid the misleading conclusion that this 
RV change is a disturbance by a companion. A feature of the 
secular acceleration is that this effect is only geometrical, not 
physical, and always positive (v, > 0). With the high astrometric 
precision of Hipparcos, this effect can be predicted with high 
accuracy for nearby stars. 

Table [3] lists the prediction for secular acceleration for all 
of our sample stars as derived from Eq. (|2]i. Some stars (e.g. 
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Table 3. Proper motion /z, parallax n (from Hipparcos, Ivan Leeuwenll2007h and secular acceleration (S A) v^if) for the M dwarfs 
sample, with the highest SA stars in the lower part of the table for comparison (not in our sample). 



Star 






n 


Vr(0 




[mas/yr] 


[mas/yr] 


[mas] 


[m/s/yr] 


Barnard 


-798.58±1.72 


10328.12+1.22 


548.31 + 1.51 


4.497+0.012 


GJ 1 


5634.68±0.86 


-2337.71+0.71 


230.42+0.90 


3.711+0.015 


GJ27.1 


485.57+2.81 


-223.13+2.11 


41.69+2.80 


0.157+0.011 


GJ 118 


978.61±2.27 


633.21+2.64 


85.87+1.99 


0.364+0.008 


GJ 160.2 


51.90±1.24 


-780.04+1.57 


43.25+1.61 


0.325+0.012 


GJ 173 


-225.75±1.94 


-192.57+1.88 


90.10+1.74 


0.022+0.000 


GJ 180 


408.07±2.49 


-642.82+2.06 


82.52+2.40 


0.161+0.005 


GJ 190 


502.99±1.32 


-1399.76+1.52 


107.85+2.10 


0.471+0.009 


GJ218 


768.13±1.39 


-123.17+1.54 


66.54+1.43 


0.209+0.005 


GJ229 


-137.09±0.50 


-713.66+0.81 


173.81+0.99 


0.070+0.000 


GJ263 


-123.73±3.68 


-814.92+2.87 


62.41+3.16 


0.250+0.013 


GJ357 


136.67±1.53 


-989.13+1.41 


110.82+1.92 


0.207+0.004 


GJ377 


-1096.84±2.11 


647.29+2.63 


61.39+2.55 


0.607+0.025 


GJ422 


-2466.98±2.87 


1180.09+2.17 


78.91+2.60 


2.178+0.072 


GJ433 


-72.5 1± 1.49 


-851.92+0.88 


112.58+1.44 


0.149+0.002 


GJ477 


-101.39±2.28 


-697.57+1.75 


52.67+3.05 


0.217+0.013 


GJ510 


-452.04±2.55 


-104.79+1.55 


59.72+2.43 


0.083+0.003 


GJ620 


-348.40±2.35 


-675.73+1.87 


60.83+2.06 


0.218+0.007 


GJ637 


-480.59±1.17 


-529.12+1.86 


62.97+1.99 


0.186+0.006 


GJ682 


-708.98±2.55 


-937.40+1.88 


196.90+2.15 


0.161+0.002 


GJ739 


153.46±2.91 


-495.17+1.98 


70.95+2.56 


0.087+0.003 


GJ817 


-918.58±3.68 


-2038.13+2.50 


52.16+2.92 


2.202+0.123 


GJ821 


713.47±2.82 


-1994.64+0.95 


82.18+2.17 


1.255+0.033 


GJ842 


888.07±1.31 


-125.56+1.16 


83.43+1.77 


0.222+0.005 


GJ855 


587.72±2.17 


-376.75+1.32 


52.22+2.17 


0.214+0.009 


GJ891 


717.43±3.40 


22.12+2.69 


62.17+3.27 


0.190+0.010 


GJ911 


-42.36±3.32 


128.87+2.44 


41.22+2.64 


0.010+0.001 


GJ 1009 


62.95±2.45 


-196.98+2.46 


55.62+2.32 


0.018+0.001 


GJ 1046 


1395.67±1.72 


547.16+2.52 


71.06+3.23 


0.727+0.033 


GJ 1100 


118.29±2.67 


-498.47+1.30 


34.57+2.79 


0.174+0.014 


GJ 3020 


-44.73±4.05 


-237.17+3.70 


43.89+4.39 


0.030+0.003 


GJ 3082 


104.14+2.04 


311.53+1.81 


60.38+1.81 


0.041+0.001 


GJ 3098 


-589.19+1.86 


-887.83+1.24 


55.98+1.91 


0.466+0.016 


GJ3671 


-603.54+1.64 


-296.43+1.36 


56.38+2.04 


0.184+0.007 


GJ 3759 


-391.22+1.50 


-411.01+1.71 


58.94+2.40 


0.126+0.005 


GJ3916 


-332.19+2.90 


-352.83+2.62 


66.21+3.18 


0.082+0.004 


GJ 3973 


-9.48+2.45 


-221.00+1.72 


54.86+2.18 


0.020+0.001 


GJ4106 


31.55+4.63 


-104.82+3.59 


9.05+3.70 


0.030+0.012 
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0.030+0.002 
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0.027+0.002 
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771.64+2.60 
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Barnard's star, GJ 1, and Proxima Cen) have such high proper 
motion that the secular acceleration can be measured (4.5, 3.7, 
and 0.4 m/s/yr, respectively). This effect is depicted with a 
dashed line for Proxima Cen and for all stars with Vi > 1 m/s/yr 
in Fig.[T] The secular acceleration was subtracted before subse- 
quent analysis of the RV data. 

We note that, even though Barnard's star has the highest (an- 
gular) proper motion jd and the dependence is quadratic on /i, 
it is not the star with the highest secular acceleration. When in- 
specting high proper motion stars from the Hipparcos catalogue 
jvan Leeuwenll2007h we found four stars that have a higher sec- 
ular acceleration because of their smaller parallax (see lower 
part of Tabled, namely: GJ 451 (Groombridge 1830, G8Vp 
+ M5.5V), GJ 9511 (K2V + K2Vfe), GJ 191 (Kapteyn's star, 
IVIIV), and GJ 9371 (sdK4). 



3.2. Tests for variability and trends 

Following the recipe outlined by E ndl et all ('2002'), we per- 
formed several statistical tests to identify variability and RV 
trends in our data. First, we asked the question for each star of 
whether the observed variability or rms cr is significantly higher 
than the mean measurement error ctrv using the F-Test (and 
F - =T as F-valuefl 

O-RV 

The calculated probabilities P(F) from the F-Test are listed 
in Table |4]for each star A low value of P{F) (e.g. < 0.01, i.e. 
99% confidence, in bold face in Table HJ indicates that the ob- 
served scatter can probably not be explained with the measure- 



- We use the one-tailed F-test because we are not interested in cases 
of error overestimation. 
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Table 4. Tests for variability and trends. Listed are the number of measurements A^, time baseline T, RV scatter cr (rms), mean RV 
measurement error cr^v, and results of statistical tests. Low values of the F-test probability P{F), as well as ^'Ot'const)' ^-S- <0.01, 
mean that a constant model is improbable, hence indicate variability (values in bold fonts). In contrast, a high P{xl\ope) indicates 
that a trend is an acceptable model while a small P(Fsiope) indicates a significant fit improvement. 
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ment errors and that there is an excess variability or a trend. 
This is the case for the stars with a high rms (GJ 477, GJ 1046, 
GJ 3020, and GJ 3916), which seem to have a companion, 
probably in the brown dwarf or low-mass star regime (see also 
Fig. |2|i. GJ 190 and GJ 263 also have a high sample variance. 
They may have a stellar companion bright enough to contami- 
nate the spectrum. We deduce this from the large measurement 
error that would occur for a double-lined spectroscopic binary 
(SB2) because our data modelling is only designed for a single- 
lined spectroscopic binary (SBl). Also, because there are only 
two measurements for GJ 190 and GJ 263, they do not stand 
out in the F-statistics. Indeed GJ 263 has already been identi- 
fied as a spec troscopic binary, an d an adaptive optics image was 
presented bv iBeuzit et al.l (|2004|) . Among the stars with an rms 
smaller than 20 m/s, only GJ 229, GJ 357, and Proxima Cen have 
an rms that is significantly greater than the measurement error. 



A similar test for RV variability is to determine the good- 
ness of fit for a constant model, i.e. calculating the;t'^ above the 
weighted RV mearQ and deriving the probability from the x^- 
distribution. These probability values P(xlon!,d mostly lower 
than P{F), and therefore there would be more variable stars ac- 
cording to our P < 0.01 criterion (Table|4]i. 

As long-period planets can cause a trend in the RV, we also 
tested for this by weighted fitting of a linear slope. A high proba- 
bility of the resulting A'siope indicates that this is an acceptable fit 
(on the contrary, a low probability indicates remaining variabil- 
ity). But it is also informative to compare ;i'^j^pg with the;^'^Qj,j( of 



^ When each measurement has the same error ctrv, one gets 



2 ^ const 
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the weighted mean via the F- value dCumming et alJlT999l) 



/^slope = (A^-2)- 



2 2 
, /^constant A'; 



slope 



Table 5. Test for periodicities in the RV (GJ 190 and GJ 263 
excluded): the best period P with its and bootstrapped FAR 



y2 

^ slope 

A high Fsiope-value indicates a fit improvement whereas a low 
probability P(Fsiope) shows this improvement to be significant 
and that it is probable not due to noise. 

In Table |4] the stars with a significant trend are marked in 
bold face. Besides the companion hosting M dwarfs GJ 477, 
GJ 3020, and GJ 3916, these are the stars Barnard's star, GJ 1 18, 
GJ 229, and Proxima Cen. The significant trend of 1.4 m/s/yr 
for GJ 229 may be caused by the wide T7p brown dwarf com- 
panion GJ 229 B, as one can show with a rough estimate. 
Assuming a circular orbit and the projected separation of 45 AU 
(a = 7.7" lNakaiimaet al.l[l995L n = 0.173") as the orbital ra- 
dius to GJ 229, the orbital period is of the order of P = 400 yr. 
Then a 30Mjup companion can cause a velocity amplitude of 
K - 170 m/s (Eq. ([T]i) and the maximum RV change is RV - 
Ky- - 2.7 m/s/yr, about twice as large as the observed trend. 
However, the low P(xliope^ indicates there is still some variabil- 
ity after trend subtraction. 

3.3. Periodogram analysis 

To test for periodicities in the RV data, w e computed the general- 
ized Lomb-Scargle periodogram (GLS, IZechmeister & Kiirstej 
12009 ). which is the equivalent of fitting sine waves including an 
offset. The adopted period search interval ranges from 2 d to the 
time baseline T of each data set. Note that 2 d will mostly exceed 
the average Nyquist frequency. However, searching at frequen- 
cies highe r than the Nyquist frequency is possible for irregular 
sampling (lPeltll2009l) . False-alarm probabi lities (FAF) were de - 
termined by bootstrap randomization (e.g. iKiirster et alJI 19971) . 
In this method random data sets are generated from an original 
data set by shuffling the RVs while retaining the observing times. 
For each random data set, the GLS periodogram was computed 
and searched for its maximum. The FAP is then given by that 
fraction of random data sets having a periodogram power higher 
than the original one. For each star we generated 1000 random 
data sets which resolves FAP > 10"^ and is sufficient to decide 
whether the FAP is below our threshold of 0.01. 

Table |5] shows that GJ 4106 and GJ 1046 have a FAP 
marginally lower than 0.01. The probable brown dwarf compan- 
ion to GJ 1046 with a minimum mass of 26.9 Mjup in an eccentric 
or bit (e = 0.28) with a 1 68.8 d period has been already published 
by iKtirster et al] (l2008h based on this UVES data set. For th e 
365 d period in Proxima Cen we refer to Endl & Kiirstej (l2008h . 
who recently analysed the RV data and identified this period as 
a 1-year alias. 

Based on data from the first 21/2 yr. IKiirster et alJ (l2003h de- 
termined two RV periods of 32 d and 45 d for Barnard's star with 
a FAP of 0.56% and 1.05%, respectively. Now, in the enlarged 
data set with a 6 1/2 yr time baseline, the second period (45 d; 
with an amplitude of 2.9 m/s) has the highest periodogram peak, 
and its FAP is now less than 0.1%, i.e. more significant (see 
Fig. m second panel for the periodogram and Fig. |5]for the RVs 
phased to this period). However, stellar activity is the probable 
cause of this variabihty and will be discussed in Sect.[ 
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Fig. 5. Radial velocity time series for Barnard's star phased to 
the 44.9 d period and the best-fitting sinusoid. 



3.4. Upper detection limits 

A Jupiter-mass planet in a circular orbit with a radius of 1 AU 
around a O.4M0-M dwarf would cause an RV amplitude of 



45 m/s (even higher for closer orbits; Eq. ([T])) and would be eas- 
ily detectable with our precision of typically a few m/s. However, 
most of our sample stars show low RV-variations and no indica- 
tion of a planet. To determine which planets in circular orbits 
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Period [d] 

Fig. 4. GLS periodograms for the Ha-index (top panel; see Sect. 13.51 1. the RV data (second panel), the RV data corrected for a Ha 
correlation (third panel), and window function (bottom) of Barnard's star 



can be excluded, we calculated detection limits in the following 
way. 

We considered the data as noise and simulated planetary 
signals by adding sine waves to the data with a range of trial 
frequencies and for 12 different equidistant phases. The time 
sampling remained unaffected. For the simulated data the gen- 
eralised Lomb-Scargle periodogram power was calculated at the 
trial sine wave frequency where the peak of the signal is ex- 
pected. If this power was below a power threshold (which cor- 
responds to an FAP of 0.01, see note below), we increased the 
amplitude of the sine wave. The simulated planet is considered 
as detected, if the power is in all 12 phases (equal to or) higher 
than the threshold. The corresponding amplitude is considered as 
the 99% detection limit. In Appendix IbI we provide an analytic 
solution to calculate this amplitude for a given power thresh- 
old and phase. This procedure was performed for a number of 
trial frequencies resulting in amplitude limits which finally were 
translated into mass upper limits (Fig. |6]l using Eq. ([TJ and the 
stellar masses from Table |2] 

It should be noted that the periodogram power p is a mea- 
sure/quantity for the FAP. By applying bootstrap randomisation 
to the original data we get a FAP vs. power relation p{FAP). No 
assumptions are necessary about the number of independent fre- 
quencies. The power value for FAP=0.01 was also taken as the 
thr eshold for t he simulated data. This modification to the method 
by lEndl et al.i (.2002) bypasses bootstrapping the simulated data 
again and is therefore more efficient and allows a dense sampling 
of the frequency. 

Here we assume that the p(FAP) relation does not change 
much when adding a sine wave, because the time sampling and 
number of measurements does not change. The only thing that 
can happen is that the rms (or above the mean) of the simu- 
lated data changes (increases) by a certain factor. But as the nor- 
malized power is invariant when the measurements are scaled 
the resulting effect is small. A compa rison of the detections lim- 
its given in Fig. |6]with those from iEndl & KiirsteJ (|2008|) for 



Proxima Cen shows quite similar results. Figure |6] also shows 
the steep increase that generally appears fo r periods longer than 
the time base (see also Gumming 2004; N elson & An gel 1998J. 

Instead of again searching the whole frequency range, only 
the power at the original sine frequency was calculated, since 
one can expect to find the simulated signal there. This saves com- 
putational effort and is more conservative, because we exclude 
spurious detections that exceed the power threshold with a lower 
amplitude at alias or noise frequencies. 

Figure |6] shows the results of our detection limit calculation. 
It can be seen that for Barnard's star or Proxima Cen, i.e. stars 
with low masses and many measurements, the detection limit 
reaches down to a few Earth masses for close-in circular or- 
bits and even within their habitable zones (HZ). Both stars have 
a priori frequencies with a FAP<0.01. These frequencies were 
excluded from the detection limit calculation (Barnard's star: 
36.1 - 36.4 d and 44.6 - 45.1 d; Proxima Cen: 295.0 - 313.11 d 
and 347.1 - 392.5 d). The HZ is de picted for each star derived 
fromFig. 1 5 of Kasting etalj (ll993t), whereas we used the stellar 
masses from Table|2] For an M dwarf with a mass of O.3M0, the 
HZ is beyond 0.1 AU. For several stars we can exclude planets 
with a few ten Earth masses in their HZ and Jupiter-mass planets 
(IMjup = 318MEanh) up to a few AU. 

3.5. Correlation between RV and Ha index? 

Stellar activity can affect the measured RV. The Ho; line is an 
indicator of stellar activity (the only available one in the UVES 
spectra). In Proxima Cen as an active flare star, the Ha line is 
in emission in contrast to low-activity Barnard's star However, 
the Ha line is variable in both cases. Therefore investigating the 
correlation between RV and variations in the Ho- line may be 
useful for correcting RV for stellar activity. 

Kiirsteretal. (2003) report a correlation of this type for 
Barnard's star. We checked this again with the now available 
longer data set for Barnard's star, as well as for the other stars. 
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Fig. 6. Detection limits for stars with at least 8 RV measurements. Pl anets in circul ar orbits above the limit would be detected with 
a probability >99%. The dashed lines show the habitable zone (HZ) (iKasting et alJfT993.) . 
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Fig. 7. Time series of the Ha-index for Barnard's star. 



As a measure of the variabi hty of the Ha line we adopt the defi- 
nition of the Hof-index bv iKiirster et al.l (l2003l) : 



Q.5(Fi + F2) 



(3) 



where Fq is the mean flux in the range of [-15.5km/s,+ 15.5km/s] 
around the Ha line (A = 656.28 nm) and Fi and F2 are the 
mean flux two reference bandpasses ([-700km/s,-300km/s] and 
[600 km/s, 1000 km/s], respectively) used for normalization. 
This index is a kind of filling-in of the Ha line (i.e for emis- 
sion 7^1) and is related to equ ivalent width (see Appendix lAb. 
Following iKiirster et alJ (l2003h . we also computed a Cal in- 
dex [-H441.5 km/s, 472.5 km/s] for comparison. The Cal Line 
(A - 6572.795 A) is expected to be stable. 

Figure I2] shows the variation in the Ha index for Barnard's 
star with time. In Fig. |8] the RV is plotted agains t the Ha in- 
dex. T he flare event has already been described in iKiirster et al.l 
(l2003h . does not seem to have any effect on the RV, and is ex- 
cluded from the correlation analysis. No further flare was de- 
tected in the longer data set. We calculated a correlation coeffi- 
cient of r = -0.42 for the whole data set. Because the new data 
(JD>2452600) have a less pronounced anti-correlation with a 
coiTelation coefficient of r = -0.25, the coii'elati on decreases 
compared to the old data set where r - -0.50 ( Kiirster et al.l 
12003 ). but the anti-correlation is still present and becomes more 
significant due to the longer baseline. 

We performed a period search (Fig. |4] top panel) for the 
Ha index of Barnard's star (again flare event is excluded). With 
the new data, a lOOOd period dominates the GLS periodogram 
(upper panel). There is also some power at 44.5 d that was 
found to be the dominant period in the RV data (second panel). 
Therefore it is probable that the 45-day RV period is caused 
by stellar activity rather than by a planet (third panel). When 
we subtract the correlation from the RV data, the rms reduces 
slightly from 3.35 m/s to 3.09 m/s. The FAP of this period for 
the coiTected RV data is only 5.5% compared to < 0.01% of the 
uncoiTected RV data (Table |5]l. A similar analysis for the other 
starfl only yields a significant correlation for GJ 433, GJ 821, 
and GJ 855 (Figs.lMTTTi. AU stars show no significant RV-Cal in- 
dex correlation. 

Finally, we investigated the question of whether more active 
stars show more RV excess scatter. Figure [12] compares the RV 
scatter (rms from Table |4]l and the relative Ha line index scat- 
ter for all M dwarfs from our sample and demonstrates that no 
correlation can be seen. 
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Fig. 8. Correlation of the Ha-index with the RV for Barnard's 
star (dots with error bars, - 0.42). For comparison the Cal- 
index (crosses). 
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Fig. 9. Correlation of the Ha- and Cal-index with the RV for 
GJ 433 (rH„ = -0.40). 
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The 6 six stars with companions were excluded. 



Fig. 10. Correlation of the Ha- and Cal-index with the RV for 
GJ 821 {rna = -0.49). 



4. Discussion 

We browsed the available literature and catalogues for known 
binarity of our M dwarfs, which may have an impact on 
planet for mation. Well known are the subst ellar companions t o 
GJ 229 (Nakaiim a etaDll995h and GJ 263 (iBeuzit et al.ll2004l) . 
Proxima Cen constitutes a widely separated common proper mo- 
tion pair with a Cen A-t-B. 

For GJ 477 and GJ 3916, binarity is indicated in the 
Hipparcos catalogue. GJ 477 has the double/multiple systems 
flag "X" (stochastic solution, probable astrometric binaries with 
a short period). GJ 3916 is listed with the flag "G" (acceleration 
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Fig. 11. Correlation of the Ha- and Cal-index with the RV for 
GJ 855 (rii„ = 0.62). 
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Fig. 12. Comparison of the RV scatter (rms) and Ha variability 
for all M dwarfs with an rms < 20m/s. No correlation can be 
seen. The very active star Proxima Cen shows only very low RV 
excess. 



or higher order terms). The large amplitudes seen in the RVs for 
both stars confirm this. 

Other stars with an entry are: GJ 620 ("G"), GJ 4106 ("G"), 
and GJ 433 ("O", i.e. orbital solution). For these stars we found 
only low RV variation. This is at least controversial f or GJ 433 
(HIP 56528). The announced period is 500 d (see also lBernsteinI 
119971) and the amplitude is expected to be several hundred m/sQ 
At 1 AU our detections limit reaches down to 0.2Mjup (Fig.|6]l. 

To our knowledge the rest of our sample so far has no dis- 
covered companions. Companions are ex plicitly exclud e d with 
near-infrared speckle interferometry by iLeinert et al.l (Il997h 
for GJ 1 and GJ 682 {AK = -4.5 mag at 1-10 AU), as 
well as GJ 891 with infrar ed coronagraphic imaging by 
[McCarthy & Z uckerman (2004) (>30Mjup at 140-1200AU). 

Kiirster et al. (2003) discussed how several stellar activity 



phenomena, such as spots, plages, or convective RV shifts, might 
affect the RV measurements. A linear RV-Ha anti-correlation 
that is present in Barnard's star could be a result of a convec- 
tive redshift caused by plages that suppress blueshifted con- 



^ The parameters for the circular orbit (e = 0, a> = 0) are P = 
500 ± 32 d, To = 2448 402 ± 28 d, aq = 4.27 ± 2.04 mas, i' = 54 ± 2, 
and n = 346 ± 22deg. With the parallax n = 112.58 ± 1.44 the semi 
major axis of the photocentre is uq = 0.0379 AU. Assuming ~ ^\ 
the amplitude would be Ki = sin;^^ = 668 m/s for the circular or- 
bit. With mi = 0.48Mo the secondary mass is in2 = 19.8Mjup calcu- 

(m sin i}^ 

~ (m+M)2 



lated from the mass function /(m) 



= ^(^v^^)^ i.< 



ismi = K- ^{m + Mf- 



p 



vective flows. We find such an anticorrelation for GJ 433 and 
GJ 821, while GJ 855 exhibits a positive correlation (Figs.l9l-[TTTl. 
Bonfils et al. (2007) presented with GJ 674 (M2.5V) an example 
where the RV and Ha variations seem to be phase-shifted. The 
correlation thus differs from a linear one and looks like a loop. 
One would expect such behaviour for a plage rotating with the 
surface. 

Flare events, such as the one observed in Barnard's star, do 
not seem to take part in the correlation. This agrees with the 
fact that no significant and strong corr elation is observed in the 
flaring M dwarf Proxima Cen (see also lEndl & Kiirster 2008.) . 



5. Conclusion 

Within the sensitivity provided by our RV precision of a few 
m/s we have not detected any planets around our sample stars. 
Most of the M dwarfs exhibit only low RV variations and some 
of them have a measurable secular acceleration due to their high 
proper motion. 

We have discussed two effects on the RV that are not caused 
by planetary companions. First, the secular acceleration, is a per- 
spective effect and causes always a positive RV trend. This ef- 
fect can be corrected easily. Vise versa, this can be seen as an 

independent measurement of the astrometric quantity ^, i.e. a 
confirmation for the ratio of squared proper motion and paral- 
lax (the absolute RV is not required). The second effect are RV 
variations caused by stellar activity. This is likely the case for 
the 45 d period we found in the RV data of Barnard's star The 
RVs correlate with the Ha index. Such a correlation was found 
here only for a few M dwarfs and therefore no conclusions can 
be drawn in general. 

As a by-product of our survey we have identified 6 M dwarfs 
with low-mass companions, four of them (GJ 477, GJ 1046 and 
GJ 3020, and GJ 3916) are brown dwarf or low-mass stellar can- 
didates and two are spectroscopic binaries (SB2: GJ 263 and 
GJ 190). Follow-up RV observations will yield the orbital pa- 
rameters and the lower limits for the companion masses m sin /. 

Our detection limits demonstrate that we can exclude gi- 
ant planets with IMjup up to 1 AU for half of our M dwarfs 
and that we are in principle capable of discovering planets with 
a few Earth masses in the habitable zones of M dwarfs with 
VLT+UVES. For this purpose an adequate number of mea- 
surements is needed to find low amplitudes of the order of the 
achieved high RV precision. The low frequency of Jupiter-mass 
planets around M dwarfs requires a large sample. The given de- 
tection limits are based on the search for single planets in cir- 
cular orbits. These limits would be higher for eccentric orbits or 
multi-planet systems. 

Even if planet detections are more spectacular, it is also 
important to report non-detections, which are required to esti- 
mate the planet frequency. Our non-detections of planets sup- 
port the increasing observational evidence of a lowe r frequency 
of Jupiter-mass planets around M dwarfs. Endl et al.l (|2006) esti- 
mated a frequency of » 1 % or less up to 1 AU orbital radius for 
Jupiter-mass planets around M dwarfs compared to 2.5% for so- 
lar like stars. This comparison cannot be done yet for low-mass 
planets because they are much harder to detect (and even more 
so for G stars), which introduces observational biases. However 
low-mass planets s eem to be quite frequent around M dwarfs 
dBonfils et al.ll2007l) . 

Expansion of current M dwarf planet searches will allow 
more precise determination of the true frequency of giant plan- 
etary companions to this type of stars and lower the detection 
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limits. Photometr ic surveys for tran sits like the recently started 
project MEarth ( llrwin et aTl l2009l) monitoring 2000 nearby 
M dwarfs and microlensing projects will also c ontribute. The 
disco veries of plarietary sy stern s around GJ 876 jPelfosse et al.l 
Il998t iMarcv et al.lll998l I2001I iRivera et al.ll2005h and GJ 581 



( Udrv et al. 20071) show that planets do exist around M dwarfs. 
This promises further discoveries of low-mass planets in the fu- 
ture with high-precision RV surveys. 
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Appendix A: Relation between index and equivalent 
width 



The equivalent width is defined as 



EW 



Fq - Fi 



where 6Ai is the width of the i-th pixel in wavelength, and Fi and 
are the flux and continuum flux in the i-th pixel, respectively. 
The equivalent width EW is measured in terms of wavelengths. 
In a normalized spectrum the continuum flux is constant = 
Fc resulting in 



EW^AA 



El 

Fc 



where A/l = 2 ^'^i is the considered wavelength range. 

The index / as defined (for Ha) in Eq. (O is, on the other 
hand, a dimensionless measure. It is normalized by the mean 
flux F - -^Yii dA/Fj taken from reference ranges instead of the 
continuum, which is sometimes diflicult to estimate, in particular 
for M dwarfs with their ubiquitous absorption lines. 

In the case that the reference regions are estimated as contin- 
uum, i.e. F\ - F2 - Fc, the index becomes t = j^jj Yii^AiFi 
and is related to EW as 



1 



EW 

~Ka' 



Note that an absorption line is indicated by EW > (0 < / < 1) 
and an emission line by EW < (/ > 1). 

Appendix B: Response of the GLS periodogram 
when adding a sine wave 

The definition of the generalized Lomb-Scargle peri- 
odogram (GLS) is (we u se the notation introduced in 
IZechmeister & Kursteiji2009l) 



1 SS YC^ + CC-YS^-lCS -YC YS 



ce ss -CS^ 



whereas abbreviations are weighted covariances for data y,, 
sine, and cosine terms (e.g. YC = 2 ^lyt cos wf; - 2 ^tyi ' 
2 Wi cos Lotj = 2 Wj{yi - y) cos wf,-; w/ are normalised weights). 

When introducing a parameter r defined by tan 2t = cc'-ss 
and replacing f,- b y t,- - f ,■ - r (resulting in CSr - 0) the GLS can 



U p y Tj 

be written as (see lZechmeister & Ktirste 



■stetl (12009b 



for details) 



YCi 



YS' 



CCr SSr 



in a form very similar to classical Lomb-Scargle periodogram. 
While the first formulation can save some computational effort, 
the use of the second is more elegant for our purpose. 

By adding a sine wave with frequency wq to the data y,, we 
generate new data x,- = y,- + acoscjoti + fesinwof; or x,- - y, + 
Qj cos (i)QTi+bj sin UQTj. This results in a new periodogram Px{<^) 
for the new data set x 



xcl xsl 

+ 



CCr SSr 
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Because the times are not changed, this neither affects the pa- 
rameter T or the sums that only depend on the time sampling 
(CCt and SSr or CC, SS, CS , and D, respectively). When 
adding the sine wave the mean changes 

X — ^ WiXi — ^ Wiji + a ^ Wi cos a>oti + b ^ w,- sin LOQti 
= y + aCo + bSi). 
The sums YC and YS change as follows 

XC — ^ Wi(xi — x) cos a>ti — YC + aCCo + bCSo 
XS - ^ Wi(xi - 1) sinwf/ -YS+ oCqS + bSSo 

= ^ Wiiyi -y + a(cos wof,- - Co) + b(sm wo?/ - So)f 

^YY + fl^CoCo + b^SoSo + IgYCq + IbYoSo + labCoSo- 

(B.l) 

As described in Sect. 13.41 it is more conservative to scan the 
power only at <y = wo instead of the whole frequency range. 
This bypasses the calculation of sums comprising two differ- 
ent frequencies and is an enormous simplification. The sums of 
simulated data can be expressed by the sums for the original 
data. Therefore is not necessary to repeat the summation for the 
simulated data. When using the notation with t this becomes 
XCr ^ YCr +arCCr and XSr^YSr + OrSSr bccausc CSr^O- 
The power response at the frequency where the sine wave 
was added is 



YC^ YS ^ 

— f + a^CCr + 2aYCr + +b^SSr + 2bYSr 

CCr SSr 



1 r 1 (l-py)YY 

= — \pyYY + XX-YY\ = 1 - — — . 

XX ^' J XX 

Now we are interested in the amplitude A or variance XX that 
is required for a given phase if to produce a desired power thresh- 
old PxioS), which corresponds to an FAR The required variance 

is XX = YYp^. 

We obtain the required amplitude by solving the quadratic 
equation resulting from Eq. ( IB. lb {a - A cos ip,b - A sin ip) 

= a^CC + b^SS + 2aYC + 2bYS + 2abCS - (XX - YY) 
- A^(CC cos^ (fi + SS sin^ (p + 2CS cos ip sin (p) 

+ 2A{YC cos ip + YS sin ip) - {XX - YY) 
= qAl- +2l3A-y 

with the substitutions a - CC cos^ ip + SSsirP'ip + 
2CS cos tp sin ip, jS - YC cos ip + YS sin ip, and y - XX — YY = 
YyE±iPl_ This leads us to the amplitude 

A{(jj,ip) - -— + 

a {-) 

which we consider as the amplitude detection limit for a fixed 
phase and for a given power threshold and which can be 
expressed by GLS sums for the original data. Probing a set of 
phases ip, we finally choose max^ A{(jj, ip). 

The second of the two solution is rejected, because we de- 
mand positive amplitudes (A > 0, ^ e [0, 360°). The terms 
a and y are always positive: y - YY ^'^'^- > (as long as 

Px — Py > 0) and a - CC cos^ ip + S S sin^ ip + 2CS cos ipsimp - 
CCrCOS^ip + SSrSin^ip > 0. 




